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I. INTRODUCTION 

Polarized Drell-Yan processes are a promising ground for the study of nucleon structure m, |2| . If the transverse 
momentum of the lepton pair is detected, polarized Drell-Yan processes can be used to probe [3, |jj various transverse 
momentum dependent (TMD) parton distributions, which have received considerable interest recently (e.g., see 0- 
H[ and reference therein). Compared with semi- inclusive dee p in elastic scattering (SIDIS), which has been used 
intensively to study TMD distributions in the last decade j9l-t2Cj|. Drell-Yan reactions have the feature that only 
parton distributions are involved, that is, there is no hadron detected in the final state. 

Theoretically, the TMD factorization for Drell-Yan process at low transverse momentum has been established, and 
also the complete leading- twist Drell-Yan structure functions for spin- 1/2 hadrons beams have been given in Ref. [2l| . 
Previous studies on Drell-Yan processes based on TMD framework mainly focused on the contribution at leading 
power of 1/Q, where Q is the invariant mass of the lepton pair. This motivate us to present a full expression for the 
Drell-Yan process at twist-three level with dilepton transverse momentum kept unintegrated, which is the main goal 
of this paper. We note that at order 1/Q, results contributed by the (transverse momentum) integrated distributions 
Jt(x), 9t{x), hr,(x) and h(x) have been worked out in Refs. [22],[23j]. We will consider both polarized and unpolarized 
scattering of spin-1/2 hadron beams. We find that at order 1/Q there are sixteen transverse momentum dependent 
structure functions for the Drell-Yan process, which can be expressed as a convolution of twist-two and twist-three 
TMD distributions. 

Experimentally, a number of polarized Drell-Yan programs have been proposed at several facilities [23 - l3l| and 
some of them could be realized in the near future to provide the first polarized data. The twist-three contributions 
can be potential experimental observables and may be accessible in certain kincmatical regions. The interest on the 
twist-three contributions also comes from the fact that they are related to the quark-gluon correlation inside the 
nucleon [HI, [33|, which is still not understood yet. 

We need to emphasize that the approach in this paper is based on the assumption that the framework of TMD 
factorization is valid at order 1/Q. The same approach has been applied in Ref. [U, [3{| to calculate the complete 
leading- twist and subleading- twist observables in SIDIS (for the production of spin-0 hadrons) , where ten twist-three 
TMD structure functions have been found. Therefore our twist-three results should not be compared with the twist- 
three mechanism I36|-l 4ll in collinear factorization that has been applied to study the single-spin asymmetries in 
Drell-Yan processes [4l| . 

The remaining content of the paper is organized as follows. In Section. II we introduce the formalism needed in the 
construction of TMD twist-three observables. In Section. Ill we present the complete expressions for the differential 
cross section of Drell-Yan process with dilepton transverse momentum unintegrated at order 1/Q. We summarize the 
paper in Section. IV. 



II. FORMALISM AND KINEMATIC AL SETTINGS 



The process we study is 



h 1 (P 1 ,S 1 ) + h 2 {P 2 ,S 2 ) ->l(l) + £(l')+X. 



(1) 
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Here we consider only the electromagnetic interaction. The notations P, and Si are the four-momenta and spins of 
the hadron beams which can be decomposed as 

M 2 

P? = Ptn$ + -±n*, (2) 
"1 

= ^n1 + P 2 -nt, (3) 

S 2 = - 1 r i —n'l-——n% + S^ T , (5) 
M 2 P 2 

where n + and n_ are two light-like vectors expressed in the light-cone coordinates, in which an arbitrary four-vector 
a is written as {a^,a + ,ax}, with = (a° ± a 3 )/-\/2 and ar = (a 1 , a 2 ). Then one can define following transverse 
tensors 

gff = tf» - n \n v _ - ntn\, ^ = d a,pa n +p n-„. (6) 

We will apply the parton model to study the Drell-Yan pair production. In this model, the leading contribution is 
from the annihilation of the quark and antiquark from each proton: q(ki)q(k 2 ) — > 7* — > ££. The momenta of the 
quark, antiquark and the virtual photon can be decomposed as 

fei = x x P+ n+ + (fcl + ^ t) n_ + fciT, (7) 

X\P{ 
(1.2 , ^2 \ 

k 2 = x 2 P 2 ~ n_ + V 2 2T> n+ + k 2T , (8) 



x 2 P 2 

q = Q= n+ +Q=r£+ 9r , (9) 

where Q 2 = q 2 , and we limit our study to the region Q\ = q 2 ^ — —q 2 -, <C Q 2 , and thus the intrinsic transverse 
momenta of quarks play significant role. 

The angular distribution of the Drell-Yan cross section is usually expressed in the dilepton rest frame (Fig. [!}, 
which can be defined by introducing the following normalized vectors [23j |: 

t = q/Q, (10) 
. 2xi - 2x 2 ~ 

^(l-c^A-c^fl,, (11) 

h = q T /Q T = (q-x 1 P 1 -x 2 P 2 )/Q T , (12) 

where Pj = p — q/(2xi). The parameter c represents the degree of freedom to distribute the transverse momentum 
between Pi and P 2 . The cases c = 0, 1/2 and 1 correspond to the Gottfried- Jackson frame [421 ], the Collins-Soper 
frame [43| and the u-channel frame, respectively. Using the normalized vectors i and z one can construct the perpen- 
dicular tensors as 

= g" v - tH v + z»z v , = -e» vpcr i p z a . (13) 

As shown in Fig. [1] the azimuthal angles of the dilepton and of the transverse spin of the hadrons with respect to the 
hadron plane can defined as [44| 

cose/) = -g^h^lxu, sin <p = e^h^lxu, (14) 
cos05 = — , — , sin (Bg = - ■ -, (15) 

/ —Cj2 y/ — S 
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where lj_ = g^lv I \J —g^l^lv and = ' S v . We note that the definition in Eqs.(Q3]) and (fT5)) is the same as the 
one used in Ref. [2l| . The two light-like vectors can be expanded as linear combinations of t, z and the perpendicular 
vector h 



1 
1 

72 



ir - z» - 2(1 - c ) 



Thus the transverse tensor and the perpendicular tensor arc related by 



9t 



V = 9? + 7f + (1 - 2c) + 0(1/Q 2 ). 



(16) 
(17) 

(18) 



where the symmetrization of indices is used. One can see that the differences are of order 1/Q. As we study the 
twist-three contribution, we need to keep track of these differences. 

In the rest frame of the dilepton, one can express the differential cross section of the Drell-Yan process as 



da 



d 4 qdn 2s q 4 



(19) 



where dfl = d cos 8d(f> is the solid angle of the lepton t. The notation L^ v denotes the lepton tensor which has the 
following form [22j ]: 



V 1 " = 



cos 2 6 



sin 2 Qz»z v 



l±ll + ^9±) + sin 26 zM^ 



(20) 



Here we have ignored the lepton masses and their polarization. 
The hadronic tensor W^ v can be expressed as [45[ 



W>1V = I e - J d2klT d2fe 2T S 2 (k 1T + k 2T - q T ) Tr k 1T ) 1 ^ a {x 2l k 2T )Y 
7 Q 4 + -y v $ Aa {x, ki T )-y^ a (x, k 2T ) + 7 M 4+l a Hx, k 2T h>*$ a Aa (x, k 1T ) 
7" ^ 7 Q $ a (x, k lT )^t a Aa {x, k 1T ) + 7° 4-^t a Aa {x, k 2 Th v $ a (x, ki T ) 



1 

QV2 

1 

QV2 



+ 



-q 



(21) 



(22) 



where the factor 1/3 takes into account color average, a is the flavor index and e a denotes the charge for flavor a. The 
first line in the curly brackets in Eq. (|22[) comes from the diagram without additional gluon connecting to the soft parts: 
the gauge-invariant TMD dependent quark-quark correlation function &(x, fcir) and antiquark-antiquark correlation 
function <fr(x, k 2 x)- The second and third lines in the curly brackets in Eq. (|22[) correspond to the diagrams involving 
one gluon which connects to one of the two soft parts, represented by the quark-gluon-quark correlator $ Aa {x, fcir) or 
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the antiquark-gluon-antiquark correlator <&\ a {x, fcir), with a restricted to be the transverse index. Up to twist-three 
level, the TMD correlator $(a;, fc<r) can be parameterized as (35l. |46| 



M 



IT 



M 
2P+ 



j_ e T h 'Tp D Ta 



M 



2M 

lp k To 



M 



9 75 



e T r ) p k Ta 



M 



+ h. 



ihf 



2M 



2M 



(23) 



The distribution functions on the right hand side (r.h.s.) of (1231) depend on x and k T , except for the functions with 
subscript s, where the following shorthand notation has been used [47| 



gi a {x, k T ) = S L 9il{x, k T ) 



kx ■ St 
M 



gi T {x,k T ) 



(24) 



and so on for the other functions. The functions with subscript "1" are the leading twist distributions which have 
probability interpretations. The other sixteen functions are twist-tree distributions. The calculations for eight T-even 
twist-three distributions has been carried in the diquark spectator model [48| and in the bag model [4!| . There is also 
attempt to calculate naive- T-odd twist-three distributions [Hfj, for which the light-cone divergence emerges. 

The decomposition of antiquark correlator <fr(x, fcy) can be achieved by the replacements x\ — > X2, n + <-> — > 

—6t, and by the relations = <3>^ r ' for T = j^, 207,^75, 175 and $' r ' = — $ c ' r ' for 1, 17^75 [47[, where c denotes the 
charge conjugation operation. 

The quark-gluon-quark correlator ^ Aa (x,kT) can be decomposed as [35| 



xM 
~~2~ 



k 



(/ ± - ^) if~(fT + i~g' T ) e Tpa S T (/; 
- (h s + i e s )-fT 75 + (h + ie) + [h^ - i e^) 



■ft) 



e pa k 9 
t T K>Tp°Ta 



M 



e Tpa k T 
M 

*7t + - 



[9t 



ap 



75) 



( ap . ■ ap \ 

{9 T +ie T 75 j 



4+ 
2 



(25) 



where the index a is restricted to be transverse. The functions on the r.h.s with tilde are interaction-dependent 
twist-three distributions. They depend on x and k T > except for the functions with subscript s, which are defined as in 
Eq. (|2"4")l . The last term inside the curly brackets will not be used in our following calculation and can be omitted here. 
In the Wandzura-Wilczek approximation [5l[, the interaction-dependent distributions are set to be zero. However, 
quantitative analysis [52J on the gi structure function shows that the violation of the Wandzura-Wilczek approximation 
is sizable. 

According to the equation of motion for the quark field, relations between twist-two distributions and twist-three 
distributions can be established [47| . A full list of these relations can be found in Ref. [Hj]. Here we only quote the 
ones which will be used in our later calculations: 



xf 1 - 


= xf 1 - 4 


h, 






m 


xg' T 


= xg' T + 


M hlT > 


xgT 


= X(Jt + 


9 IT + TT "IT- 






p T m 
2M* 9lT+ M Hl ' 


xgr 


= xg T - 


xgi 


= xgi + 


m , 
9lL+ M 1L ' 


xh L 


= xh L + 


p T ± m 


xKt 


= xKt — 


Pt , m 




= xJit 4- 








Pt f ± 
JphT, 


xf T 


= xf T + 



9lT, 



(26) 
(27) 
(28) 

(29) 
(30) 
(31) 
(32) 
(33) 
(34) 
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x fr 


— x It 


+ flT, 


(35) 


xfr 


= xf T - 


+■ 2M 2 hT, 


(36) 




= xft, 




(37) 


xg 1 - 


= xg 1 - 


m , 1 


(38) 


xh 


= xh + 


M2 hl ' 


(39) 



We note that the Lorentz invariance relations [22|, |U [H, H3| for integrated twist-tree distributions and their viola- 
tions I55T157I have been discussed in literature. 



III. EXPRESSIONS FOR TWIST-TREE STRUCTURE FUNCTIONS OF DRELL-YAN PROCESS 



By substituting Eqs. ([23| . (|25j) into Eq. (|22]l. and then contracting the lepton tensor and the hadronic tensor, one 
can get the expression of the cross section for the Drell-Yan process up to order 1/Q. As shown in Ref. [2l|, there are 
forty-eight independent structure functions for Drcll-Yan processes from the collision of two polarized spin-1/2 hadrons 
beams, and twenty-four of them arc leading-twist obscrvables. Here we restrict ourself to consider the contributions 
at order 1/Q. We found that the angular distribution of the differential cross section of polarized Drell-Yan process 

at twist-thrpp IpvpI pan Vtp pvnrpssprl as 



j twist-3 
LLU _ "em 

dx\dx 2 d 2 qr dfl 3Q 2 



sin 2$ | cos 4> Fuu ' 



S 1L sin d> Ff£ + S 2L sin <j> F™ * + S 1L S 2L cos F™ s ' 



I -Sit | 
\S 2 t\ 



sin(ffi>i 



\ p sin(0 Sl +ii 



sin(^s 2 



)F, 



,sin(0s 2 +<; 



UT 



■ S±L \S 2 T 

■ S 2 L \Sit 
- \S\t\ \S 2 t 

■ cos(0 Sl - 0S 2 



- COS(0s 2 + 0) 

- cos((p Sl + 0) 

cos(0 Sl + <j>S 2 

cos' 
TT 



LT 



- sin(0s 2 - 4>)F, 

■■(4>s 2 +4>) 



,sin(4>s 2 
UT 



■ COS(05 2 



\F\ 



,cos(0s 2 — « 



n cos(0 s , +<P) , // 
TL +COS(<j> Sl -<p)r TL 



LT 

,cos(0si 



)F 



,COS 

TT 



,-<p)f: 

<*s 1+ * aa+ *) +cos{lj)8i+<f>Sa _ 



\ 7-t cos (0s 1 +4>s 2 -<t>) 
)r TT 



,cos(0 Sl -0 S2+ 0) + cQg( ^ _ ^ _ ^ F co.(0 Sl -0 S2 -0)j I 



(40) 

consistent with the results 



We obtain sixteen twist-three structure functions, the angular dependences of which are 
given in Eq. (57) of Ref. [p. 

To shorten the notation we will use following combinations, since they always appear in the same way: 

' = xi ((l-c)/ + c/), J = x 2 (cj+(l-c)f 



(41) 



where / and / stand for the twist-three quark distributions, and / and / for the antiquark distributions, respectively. 
Using the EOM relations in Eqs. (|26|) to (|39l) . the structure functions for polarized Drell-Yan processes at twist-three 
thus are expressed as 

(h ■ k 1T ) (f^ ^rhit) (h ■ k 2T ) (hJ 1 - ^-hhi 



F, 



iCOS ( 

UU 



f; 



isin <j) 



UL 



-3 e 



771SIII <p 

b LU = 



ic 



(h ■ k 1T ) (g 1 - Vil + ]jf >*i h L \ - (h ■ k 2T ) (fCfi + ^} h i 



pSin(0s 2 -0) _ J_ r 
r UT 



(h . k 1T ) (jtfi + ^ hi L t\ - (h ■ k 2T ) L^f- + ^-h L hf 



(42) 
(43) 
(44) 

(45) 



f: 



sin(</>s 2 +4>) 



UT 



F 



TV 



F 



b Sl -4>) 



TV 



LL 



F\ 



LT 



F 



LT 



<<p2 + <t>) 



F 



COS(0S! -<t>) 



TL 



F 



COs(0 Sl +0) 



TL 



TT 



F 



cos(0 Q +(ji(,+0) 



TT 



2 ft-fc 



IT 



k 2 
"-2T 



hfr , M 1 hh^ j 



IT 



Mo 



Mi 



^2ft ■ k\xh ■ k2T — k\T ■ fc 2 T 



f ± flT , 9 X 9lT , /ll"^T , ^l"^ 



Mo 



M 2 



Mi 



Mi 



■5" 



2Mi /t/! + 2M 2 hih + (feiT ■ fe 2T ) 1 J1T 



/it/ 1 " 51T3 1 ^1" . h^h^ 



- 2 h-fei T -fc 



"IT 



2ft ■ kyrh • k 2 T — &it ■ &2T 



Mi 

Mf / 

giTg 1 
Mi 



Mi 



M 2 



M 2 



Mi 



A-t/ 1 
Mi 



M, 



M 2 



3 e 



(ft • fc 1T ) ( fflL + -Ih^hL) - {h ■ k 2T ) (g XL gl + ^-h L hf L 



Mi 



M 2 



Q 



2M 2 3lL 3 T + 2M 1 h L h x - (fcir • k 2T ) 



( dhhT + h iL h T TtfvT _ Kl H T 



Mf 



2ft ■ kirh ■ k 2 T — feiT ■ &2T 



^ M 2 Mi M 2 Mi 



Q 



-2Mi5T3 1L - 2M 2 ftifci + (fciT • k 2T ) 



- 2 (h-feiT -fe 



"IT 



^2ft • kixh ■ k 2 T — k\T • &2T 



V Mi 

?t ffiL , M 2 h{ T h L 
Mi Mf 

| Sit si , ftr^ii 



/ fflTgL ^T fe^L 

I Mi Mo 



Mi 

flT ft 

Mi 



^T ^1L 
M 2 



Mi 



M 2 



At f_L 
Mi 



fey /ij^ 

Mo 



(ft-fc 



IT J 



M 2 /fV /t M 2 giT ff T 



Mi 



Mi 



/ly/ll — ft, T fei 



(h.fcar) f M2Mrr- M ^{^T +hi % T + hi %± 



M 2 M 2 
4(ft • fc 2 r)(ft ■ fcir) 2 - 2 (ft. ■ k 1T )(k 1T ■ k 2T ) - (ft ■ k 2T )k\ T 

It /it St Sit ^it 



\ 2M X M 2 2M X M 2 1M\ 2M 2 I 

(4(ft • fei T )(ft ■ k 2T f - 2(ft ■ fc 2T )(fc 1T ■ fc 2T ) - (ft ■ k 1T )k 

( /it /t Sit St ^t h 1T h 



2M X M 2 2MiM 2 2M| 



T "IT 

2M| 



7 



TT 



s(<t>a-<t>b + <t>) 



(h ■ k 



IT) 



F, 



cos(0 o -0i,-0) 



TT 



■C 



M 2 frr fr _ m i9\t9t 
M x M x 

2(h ■ k 1T ){k 1T ■ k 2T ) - {h ■ k 2T )kl T ^ 

frfjr , 9t 9 it , ^t _ h^rphj, 
l 2M 1 M 2 2M X M 2 2M\ 2M\ 

Mif T fiT M{g T g 1T 



hrhi — h T hi 



(hk 



2T) 



Mo 



Mo 



hi hx — h\h T 



— ( 2{h ■ k 2T )(k 1T ■ k 2T ) - (h ■ ki T )k 2T 



/it /t i 9it9t , hrh^ T h^h^- 



<>T n lT 



\2MiM 2 2M X M 2 2M| 2M| 
In the above equations we have used the notation: 

C[w(k 1T ,k 2T )fg] = J2 e l f d2k ^ d 2 k 2T 
q J 

x <5 (2) (q T - fei T - fc 2 r) w(k iT , k 2T ) 
x [/'(si, tfr^te.A&O 
+/ ? (zi,fc 1 V)s ? (z2,4r)] . 
Also, we have applied the following combinations for certain distribution functions: 



fx(x, k T ) = fx(x, k T ) 
<7 T (x, k T ) = g' T (x, k T ) 

9t( x i ^t) = 9t( x > ^t) — 
hi(x, kj.) = hrr{x, k^) - 



k 2 



2M 2 



It { x i ^t)j 



k 2 2 
9t \ x i ^t)i 



2M 2 
k 2 



2M 2 



2M 2 



9t ( x i ^t)j 



k T 2 



(56) 
(57) 

(58) 



(59) 



(60) 
(61) 
(62) 
(63) 
(64) 



Equations (|42p to (|58| represent a main result of this work. They are complementary to the leading-twist structure 
functions of Drell-Yan processes. We note that, after being integrated over qx, the expression in Eq. (|40"]) reduces to 
the following result [H[ 



da 



twist-3 



01 em sin26»4{|S'iT|sin 



dxidx 2 dQ 3Q 2 Q 

+ \S 2T \ sin(^s 2 - 4 

+ Sil\S 2T \ cos(cj>s 2 

- S , 2l|<S'it|cos(^ Si 



-0) (Mif T fi+M., 
M 2 fi% + Mihh^j 
4>) (aI 2 gig T + Mih L hi 
Mi9T9i+M 2 hih L ^ }. 



(65) 



where the distributions are the transverse-momentum integrated version. The last two terms in the curly brackets 
were given in Ref. [22[, while the first two terms were given in Ref. [23|. However, as Jt(x) and h(x) vanish because of 
time reversal invariance (46| , the first two terms in Eq. (|65j) only receive contributions from the interaction-dependent 
twist-three distributions fx and h. The function Jt{x) is related [45|, [H| to the collinear twist-3 Efremov-Teryaev- 
Qiu-Stcrman (ETQS) function Tp(x,x) [36l - [38| . The role of the ETQS function in the single spin asymmetry in 
Drell-Yan were studied intensively in the literature [58U6l| . and has been revisited in Refs. [62| and [63| recently. It 
was found that different theoretical approaches lead to different results. For example, an additional factor 2 was found 



8 



in Rcf. |62j while factor 1/2 was gained in Ref. [63[, in comparison with the expression for An in Rcfs. [23l. |60L l61j. 
We note that the normalization of the asymmetry of our result (in the Collins-Soper frame) agrees with that in 
Rcfs. @, Ho|, ED], but not with that in Rcfs. [HI El]. 

In the following we present further comments on our result. 

• Certain structure functions in Eqs. (|42[) to (|58|) satisfy the following symmetric or antisymmetric property: 

p sin0 _ _ p sin0 p sinOs 2 -<t>) _ p sin(0 Sl -<P) 
r UL ~~ r LU ' r UT ~ r TU ' 

„sin(0 S2 +0) _ sin(0 Sl +0) cos(</> s 2 - (p) _ cos(<p Sl ~<t>) T? cos (</>S 2 +<f>) _ _ p cos(0 Sl +0) , , 

which agree with the general analysis given in Ref. [2l| . 

• The expressions of the twist-three structure functions depend on the choice of the dilepton rest frame, char- 
acterized by the parameter c. This feature is different from that of leading-twist structure functions, whose 
expressions of which are the same in different dilepton rest frames. 

• Among the 16 twist-three distributions, 12 of them appear in the twist-three structure functions for Drcll- 
Yan processes, including the new distributions g 1 - [34| and /jr (46[. The first one appears in both the single 
longitudinally and transverse polarized Drcll-Yan processes, while the later one appears in single and double 
transversely polarized Drcll-Yan processes. These new twist-three TMDs are quite interesting, since evidence 
of their existence will indicate the necessity of introducing the gauge-link direction in the decomposition of the 
correlator I34l. |4 6| . Thev represent new nucleon parton structure information that has not been explored before. 
Calculations [64l . |65( in a diquark model show that the contribution from g 1 - should be included in order to give 
a complete description of the asymmetries in longitudinally polarized SIDIS. 

The distributions e, ej,, ex and do not contribute here, since we have not considered the lepton polarization. 
These four distributions could be studied in SIDIS with a polarized lepton beam [l2|, [H, |66|, [67} . 

• Equation (|42|) shows that the combination f 1 - g) f 1 or ® h and so on can lead to a cos(f> asymmetry in 
unpolarized Drell-Yan processes, which is similar to the case in SIDIS f(38t |. Apart from the well-known cos2</> 
distributions, the measurements of dilepton production in unpolarized hadron-hadron collision also show a 
cos^> angular dependence [69l - l7l| . It has been shown that the QCD corrections can generate such angular 
distribution |72h77| | . Our study shows that there is an alternative mechanism for the cos <fi asymmetry in the 
unpolarized Drell-Yan process, due to the presence of twist-three TMD distributions. 

• In the region where the TMD framework is assumed to be valid, there is a suppression factor qr/Q for twist- 
three structure functions. Therefore the asymmetries arising from twist-three structure functions are supposed 
to be smaller than the leading-twist observables. However, the asymmetries are also directly determined by the 
size of the twist-three TMD distributions. It is not known whether there are positivity bounds to constrain 
twist-three TMD distributions, like the case of leading- twist distributions [78|. Sizable twist-three TMD distri- 
butions could lead to nonvanishi ng asy mmetries in Drcll-Yan process. Furthermore, the SIDIS measurements 
on fixed targets 0, [ll|, [lH, Gjl [Sli I79l - l8l| show that the twist-three asymmetries are not small. This also encour- 
ages the cor resp onding measurements of twist-three effects in Drell-Yan processes, especially for fixed-target 
experiments 0, (H-Hi, [H ■ 

• For each structure function, there are several combinations of twist-two distribution and twist-three distribution 
which can contribute, similar to the twist-three results of SIDIS process. This makes twist-three parton distri- 
butions more difficult to be probed in high energy processes than the twist-two ones. Further theoretical and 
experimental studies are needed to provide more constrains the size of different TMD twist-three distributions. 



We point out that our calculation is based on a generalization of the TMD factorization to the twist-three 
level. Therefore the correctness of our results relics on the validation of the twist-three TMD factorization. 
Unlike the twist-three collinear factorization, which has been widely applied in SIDIS and Drell-Yan, the TMD 
factorization formalism at twist-three level has not been established yet. The main challenge for twist-three 
TMD observables is that when one calculates the twist-three TMDs, there are light-cone divergences [5(3| for 
which it has not been understood how to control them at order l/Q. This does not necessarily means that 
twist-three TMD factorization cannot be developed. Further study is needed to overcome this difficulty. 
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IV. CONCLUSION 

Drcll-Yan process has been recognized as an important tool to study the structure of the nucleon. In this work, we 
have studied polarized Drcll-Yan processes from the collision of two polarized spin- 1/2 hadrons beams at order 1/Q, 
based on the framework of TMD factorization. We find that, among a total of twenty- four sublcading- twist structure 
functions in Drcll-Yan process, sixteen of them are at twist-three level and can be expressed as combinations of twist- 
two and twist-three TMD distribution functions. We give the complete expressions for these structure functions, for 
each of which there are several twist-three distributions that contribute. Based on our result, we point out that twist- 
three distributions can provide an alternative explanation for the cos <f> angular dependence in unpolarized Drell-Yan 
processes. The measurements of the asymmetries at order 1/Q in Drell-Yan process therefore can provide useful 
information on the twist-three TMD distributions and the multi-parton correlations in the nucleon. 
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